Characteristic spatial scale of vesicle pair interactions in a plane linear flow 
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We report the experimental studies on interaction of two vesicles trapped in a micro-fluidic four- 
roll mill, where a plane linear flow is realized. We found that the dynamics of a vesicle in tank- 
treading motion is significantly altered by the presence of another vesicle at separation distances up 
to 3.2-^3.7 times of the vesicle effective radius. This result is supported by measurement of a single 
vesicle back-reaction on the velocity field. Thus, the experiment provides the upper bound for the 
volume fraction (f> — 0.08 -J- 0.13 of non- interacting vesicle suspensions. 

PACS numbers: 87.16.D-, 82.70.Uv, 83.50.-v 



The coupling between rheological, macroscopic prop- 
erties of complex fluids, such as suspensions, emulsions, 
polymer solutions etc, and the microscopic dynamics of 
the deformable micro-objects immersed in them has been 
a long-standing problem in physics, chemistry, and en- 
gineering. The simpler the micro-object dynamics, the 
more impressive the progress in understanding rheology 
of the corresponding complex fluid. A prominent exam- 
ple is a colloid suspension, where the correction for the 
effective viscosity of a dilute suspension due to immersed 
solid spheres has been obtained about a century ago [lj. 
However, for such complex fluids, as suspension of vesi- 
cles, capsules or red blood cells, which obviously have 
great biological and industrial applications, the relation 
between their rheological properties and the dynamics 
of the deformable micro-objects still remain under dis- 
cussion due to their more elaborate microscopic behav- 
ior. For vesicle suspensions, the currently available the- 
ories deal only with the dilute regime, so that its 
rheology is described in the limit of a single vesicle by 
completely neglecting vesicle interactions. On the other 
hand, the results of the two currently available experi- 
ments on the dependence of the effective viscosity of the 
vesicle suspension fj on the vesicle viscosity contrast be- 
tween the inner and outer fluids A contradict each other 
in measurements taken at about the same volume frac- 
tion of vesicles of 0.11-0.12 @, Q. Though one of the 
experiments @] agrees with the theory 0, the other 
[6] shows a non-monotonic dependence of fj on A, con- 
trary to the predictions for dilute suspension J2|-|5]. The 
results of the experiments reported in Ref. [6j indicate 
that hydrodynamically-assisted vesicle interactions could 
be a cause for the observed qualitative difference with the 
predicted behavior. Indeed, both vertical velocity fluctu- 
ations and vesicle inclination angle in tank-treading (TT) 
motion in a shear flow were significantly affected by many 
surrounding vesicles, as measurements presented in Ref. 
Q showed. In this respect, the key question is whether a 
long-range hydrodynamically-assisted interaction of vesi- 
cles in flow can qualitatively modify the rheological prop- 
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Figure 1: (Color on line) Experimental setup. 



erties of the suspension and at what characteristic length 
scale does the vesicle interaction become significant to 
achieve this. To the best of our knowledge, the role of 
long-range hydrodynamic interaction was addressed the- 
oretically [8| and numerically [9] only for colloids, where 
interaction of two suspended spheres in a shear flow was 
studied. 

Here we report the studies of the correlation and the 
range of the interaction of two vesicles in a plane linear 
flow. By applying a novel and unique technique for vesi- 
cles, deviations of the vesicle inclination angle in the TT 
motion from the stationary value were used as a sensitive 
quantitative detector of the interaction strength between 
vesicles. It provides the characteristic length scale of the 
long-range hydrodynamically-assisted vesicle pair inter- 
action, which determines a lower bound of the interaction 
scale for a semi-dilute vesicle suspension. This method 




Figure 2: (Color) The effect of a nearly spherical vesicle 
on the velocity field. The 2D time-averaged velocity fields 
(over 122s on 7 vesicle orbits) are presented: : (top) reference 
field without a vesicle; (middle) in the presence of a nearly 
spherical vesicle (R ~ 17 /im); (bottom) difference field, i.e. 
subtraction of the reference field from that measured in the 
presence of the vesicle. Colors indicate velocity magnitude 
[fim ■ s" 1 ]; Vesicle trajectory is denoted in black. The plotted 
circle shows the vesicle size and mean position. Note that 
velocity differences exceeding 1 [im/s are significant. It was 
verified by comparing two reference flow fields. 



Figure 3: (Color) The effect of a deflated (non-spherical) vesi- 
cle on the velocity field. The 2D partial time-averaged velocity 
fields are presented; in the presence of the vesicle the averag- 
ing is taken over 2s segments (marked in red) along each of 6 
vesicle orbits (in black) : (top) reference field without a vesicle; 
(middle) in the presence of a deflated vesicle (R ~ 16.6/im, 
A ~ 0.39); (bottom) difference field, i.e. subtraction of the 
reference field from that measured in the presence of the vesi- 
cle. Colors indicate velocity magnitude \p,m ■ s -1 ]. The plot- 
ted ellipse shows the vesicle size and mean position during the 
time over which the averaging was taken. Note that velocity 
differences exceeding 1 /zm/s are significant. This was verified 
by comparing two reference flow fields. 



could in principle be applied to other deformable objects. 

A vesicle is a drop of fluid, encapsulated in a lipid 
bilayer membrane and suspended in either the same or 
different fluid. The impermeability and the inextensibil- 
ity of the lipid membrane dictate conservation of both 
volume and surface area of the vesicle. Vesicles undergo 
a TT motion in a planar linear flow at A = 1 and uj / s < 2 
with < A < 1 Here A = A/R 2 - 4tt is the excess 



area, R = (3V/4TT) 1 / 3 is the vesicle effective radius, V 
is the vesicle volume, A is the vesicle surface area, and 
to and s are the vorticity and strain rate, respectively 



At a fixed viscosity ratio, the inclination an gle 9 



decreases with increasing A and increasing uj/s [ll - 14 1. 

The experiments were conducted in a micro-fluidic 
four- roll mill [lol . 15|. implemented in silicone elastomer 



(Sylgard 184, Dow Corning) by soft lithography. The 
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micro-channel was kept in a vacuum chamber for > 1 
hour before filling with vesicle solution, to avoid air bub- 
bles. The key component of this device is a dynamical 
trap, which allows a long observation time compared with 
the orbit period of vesicles in the flow. The observation 
time is limited by the fluorescent lipids bleaching. The 
flow was driven by gravity and the control parameter uj / s 
was varied continuously by variation of the pressure drop 
AP across the device in the range [0-0.01] Pa (Fig. []}. 

Vesicles were prepared in water (rjout — r\in = 
ImPa ■ sec) via electro- formation |la ). using lipid 
solution, consisting of 85% DOPC lipids (Sigma) and 
15% NBD-PC fluorescent lipids (Molecular Probes), 
dissolved in 9:1 v/v chloroform/methanol [3, EH- The 
electro-formation cell consisted of two Indium Tin Oxide 
(ITO) coated glass electrodes with a 1 mm teflon spacing 
between them. Droplets of the lipid solution were spread 
on the glasses. The solvent quickly evaporated leaving 
just the dry lipids. The cells were then filled with 
de-ionized water (Merck, R> 1M =cm). An external 
AC voltage, sine wave, 3 volts peak-to-peak, 10 Hz, 
was applied to the electrodes. The electro-formation 
cells were placed on an orbital shaker, rotating at 90 
rpm for > 8 hours. New solution was prepared prior 
to each experiment [12, [l3j]. The vesicles dynamics 
were monitored using inverted fluorescence microscope 
(IMT-2, Olympus) inside a 180 x 140/im observation 
window, parallel to the flow (Fig. [TJ. The images were 
collected using a Prosilica EC1380 CCD camera, aligned 
with the shear axis, at 30 fps, with a spatial resolution 
of 0.27 (im/ pixel. The camera was synchronized with a 
mechanical chopper on the path of the excitation beam 
to reduce exposure time. Vesicles were loaded into the 
dynamical trap via an inlet channel at the mid-plane 
of the micro-channel (of 300 fim height). To isolate 
just two vesicles near the trap center, we repeatedly 
varied AP quickly (on the order of a second) generating 
short elongation flow pulses (open trajectories), which 
expelled the vesicles located far from the trap center. 
After these perturbations and setting the fixed value of 
ui/s, we allow the flow to relax during several vesicle 
relaxation times, which is of the order of 5-10 s 
before starting the measurements. Moving along closed 
trajectories imposed by the external flow field, the 
vesicles were repeatedly approaching and moving away 
from each other (Fig. [5]c and movie in [13] )■ A special 
attention was given to having the largest cross-sections 
of the vesicles in the same plane in order to reduce 3D 
interaction effects. Image processing was based on a 
binary recognition of each vesicle used to extract the 
coordinates of its center. The relative position of the 
membrane was determined in the frame of reference of 
the vesicle using intensity variations along the radial 
directions. The subsequent analysis included elliptical 
approximation of the shape and determination of the 
lengths of the main axes, the inclination angle, and the 



center of the ellipse. Following [l2j, Il3 |. both R and A 
were calculated from the main axes. 

Before presenting our main results regarding the 
interaction of two vesicles, it is insightful to examine 
direct measurements of the velocity field in the vicinity 
of a single vesicle, demonstrating the back-reaction of 
the vesicle on the flow. The velocity measurements were 
conducted by micro Particle Image Velocimetry (/iPIV). 
In this experiment an aqueous vesicles suspension 
was mixed with fluorescent particles (Fluoresbrite® 
YG Microspheres 0.5/ito from Poly sciences, Inc.) at 
2.6 x 10 7 particles ■ ^Z _1 , allowing imaging of the flow in 
the four-roll mill trap. Images were recorded at 90.9 fr/s 
rate on a computer using the Motmot Python camera 
interface package [18], controlling a GX1920 camera 
from Allied Vision Technologies. The images were then 
filtered (Laplace filter using Gaussian second derivatives 
[lj|) and processed using Gpiv [20j at interrogation 
windows of 32 x 32 pixels (corresponding to about 
8.6 x 8.6/ito 2 ) with 50% overlap. The effect of a vesicle 
with vanishingly small A, located near the stagnation 
point of the trap and moving on a small orbit in a plane 
linear flow with 10 = 0.46 s _1 and lu/s = 1.8, on the 
velocity field was measured. Figure [2] at the bottom 
presents the difference velocity field, which results from a 
subtraction of the reference field (Fig|2] at the top) from 
that measured in the presence of the vesicle (Fig. [2] at 
the middle). It allows to study the back- reaction of the 
vesicle on the flow, and to estimate the effective distance 
of the back-reaction. The most pronounced feature 
of the field difference is the alternation of the strain 
strength in the vicinity of the vesicle, which is significant 
at scale of up to about 3i? along the strain main axes. 
The compression and stretching directions of the field 
difference show sign inversion compared to those found 
in the reference flow field, where a vesicle is absent (Fig. 
[2] at the top). This can be attributed to the vesicle's 
volume and area conservation, opposing extension and 
compression. It can be also viewed as a result of the 
tank-treading motion of the membrane, which must 
move at constant tangential velocity. Indeed, it was also 
observed in a snapshot of velocity field around a vesicle 
in 2D numerical simulations [2l| . Very similar image 
of the disturbance velocity field around a solid sphere 
in a shear flow was recently obtained by numerical 
simulations [jj. 

Similar analysis of the difference velocity field around 
a deflated (non-spherical) vesicle moving on a small orbit 
in a plane linear flow with ui — 0.46 s _1 and u/s =1.4 
after partial time averaging on small part of its six 
orbits as shown in the vesicle center is presented in Fig. 
[3] The pronounce difference with Fig. [2] is symmetry 
breaking of compression and stretching directions in 
the disturbance velocity field: the stretching direction 
deviates from the strain direction by about 7r/12. 

We conducted two types of observations of vesicle 



4 



(a) 



(c) 



^ "^fk^T^ 



(b) 



'.z/sec y.zjsee lu.jbsec n.^isec u.wsec I4.iusec 



Figure 4: (Color on line) Snapshots of vesicle pair interactions, a) A ~ 0, 1, R — 13.2, 9.2/^m, uj/s = 1.33. b) A 
0.88, 0.77, R ~ 10.1, 11.9pim, ui/s — 1.52. c) Trajectories of interacting vesicles inside the dynamical trap at to/s = 1.52 
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Figure 5: (Color on line) The upper plot shows a time series 
of 59, the deviations of the vesicle inclination angle from its 
stationary value (blue squares), for a deflated vesicle along 
its orbit as it interacts with a spherical one. The green dots 
denote SO of the same deflated vesicle after removal of the 
spherical vesicle. The lower plot shows a time series of the 
distance between the centers of the interacting vesicles. 



pair interactions: interaction of spherical and deflated 
vesicles in a pair and interaction of two deflated vesicles 
in a pair. First we observed interaction of a deflated 
vesicle in the range of parameters A G [0.7, 1] and 
R G [6, ll]/im, placed near the stagnation point of 
the flow, and a vesicle with vanishingly small A and 
R G [ll,17]/im rotating around the former, all together 
nine such pairs. By choosing one spherical vesicle in 
each pair, we reduce the number of parameters and 
simplify quantification of the problem (see Fig. 2] a and 
movie in Ref. [l3|)- We took sufficiently long time series 
to observe several periods of the vesicle interaction. The 
latter is characterized by 89 = 9 — (9), the deviation 
of the inclination angle of the deflated vesicle 9 from 
its mean stationary value (9) fill . Il3| and is typically 
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Figure 6: 86 versus a/n. Two data-sets are for A ~ 0, 1, R : 
13.2,9.2 fim, lo/s = 1.33 - circles, and A ~ 0,0.82, R ■ 
12.5, 7.5 fim, tj/s = 1.27 - squares. 



accompanied by the vesicle deformation. The observed 
89 due to the vesicle interactions are up to one order of 
magnitude larger than due to thermal fluctuations. 

To relate the deviations of the vesicle inclination 
angle from its stationary value 89 for a deflated vesicle 
along its orbit, as it interacts with a spherical one, to 
the distance between vesicle centers we added two plots 
in Fig. [S] To determine 89 as a function of time for 
the same vesicle, 9(t) of the same deflated vesicle was 
measured in the presence of a spherical one and then 
after its removal (Fig. at the top). At the same 
time, the distance between the centers of the interacting 
vesicles was measured (Fig. [S]at the bottom). 

In order to relate the data to the structure of the 
strain rate of the velocity field around a vesicle (Fig. 
[2]), we plot 89 versus a, the angle between the line 
connecting both vesicle centers and the horizontal 
direction (see inset in Fig. In the two different 
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Figure 7: I R. The data are taken over many 

cycles in time series of 9 vesicle pairs. Inset: The same data 
on logarithmic scale. Dash line is the fit ag ~ (dmin/R) 13 with 
= -2.65±0.2. (see also Fig. 1 in where different vesicle 
pairs are shown by different symbols) 



data sets presented in Fig. [51 where each data set 
contains 3-4 interaction events, we observe a repetitive 
pattern, periodic in a/n. The pattern is characterized 
by maxima of 59 occurring in the stretching direction 
of the velocity field difference at a/n = —1/4 and 
3/4, while 59 minima are shifted from the compression 
direction at a/ir = —3/4 and 1/4 (see Fig. The deep 
minima of 59 result from a short-range vesicle-vesicle 
interaction, which is accompanied by strong vesicle 
deformation (see Fig. @] a) at distance d/R w 1 and 
bounded by a width of a fluid lubricating layer, where 
d/R is the normalized distance between vesicle centers, 
and R is the radius of the nearly spherical vesicle. The 
short-range interaction was studied experimentally for a 
vesicle @ and droplet 22 1 pairs and numerically for a 
capsule pair [23j|. The maxima of 59 occur as the result 
of a long-range hydrodynamically-assisted vesicle-vesicle 
interaction at d/R determined above. It has not been 
studied before for any soft objects and appears to be the 
main interest of our investigation here. 

Figure [7] presents ag , the rms of the deflated vesicle 
inclination angle 9 (inset in Fig. 5) as a function of 
d/R in the vicinity of its minimal value d m i n /R, the 
relevant parameter in the problem. As seen in the plot, 
ag reduces with d m i n /R and reaches background noise 



level of the order 0.02 -=- 0.03 rad at d min /R w 3.2 -=- 3.7 
(in Fig. 1 of [l3| the same data are shown, and different 
vesicle pairs are presented by different symbols). This 
noise level, measured independently for a single vesicle 
in the four-roll mill, is attributed to thermal fluctuations 
and non- uniformities of the flow (see Fig. [5j. The same 
data, shown in the inset in Fig. [7] on a logarithmic scale, 
is fitted by ag ~ (d min /R)P with (3 = -2.65 ± 0.2. It 
should be compared with a long range hydrodynamic 
particle-particle interaction in a shear flow, where the 




O 



-10 -5 5 101 

x[sec] 



0.5 to 
E 

O 

,0 



d . /<R> 

min 



Figure 8: Characteristic dynamics of two deflated vesicle pair, 
A = 0.876, 0.737, R = 10.3, W.lSfim. a) d versus time, b) 59 
versus time, c) C(r) versus r. d) C max (r) versus d min /{R). 



presence of a second particle changes the induced strain 
rate of a test particle as d~ 3 for unbounded and dT 2 
for bounded suspensions jH, The spatial scales, at 
which the vesicle pair interaction reduced down to the 
noise level, corresponds to a vesicle volume fraction of 
0w 0.08 -=-0.13. 

In a second, general case of two interacting deflated 
vesicles (see Fig. @] b,c) more parameters characteriz- 
ing both the vesicle and flow geometry (see Fig. [3]) 
are involved, which makes the situation less quantita- 
tively tractable. In this case a correlation of the incli- 
nation angles of the vesicles near a local minima of the 
distance between vesicle centers was observed (Fig. [5] 
a,b). We investigated eight pairs of vesicles with simi- 
lar physical parameters in the range A 6 [0.5,0.9], R 6 
[6, 12][im and found high correlation, measured via the 
peak C max {T) of the normalized cross-correlation func- 
tion C(t) = (59i(t)59 2 {t + r)) t /ag 1 ag 2 for nearly collid- 
ing vesicles (see Fig. [5]c), followed by its rapid decrease 
in the region of d min /(R) w 2-4-3.3 (Fig. g]d). Here 
(R) is the arithmetic average of the radii of the two vesi- 
cles in pair. Thus the characteristic length scale defined 
from Fig. [8]d, where C m ax(T) reaches the noise level, is 
dmin I '(R) ~ 3.3 that corresponds to a volume fraction of 
0ra 0.12 ±0.03. 

To conclude, we have found that the dynamics of a 
single vesicle can be significantly modified by the pres- 
ence of another vesicle due to long-range hydrodynamic 
interaction that is also strengthened by the direct mea- 
surements of a single vesicle back-reaction on the velocity 
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field. A direct link to the rheology of vesicle suspensions 
still remains a challenge due to the diversity of the phys- 
ical parameters and the many- vesicle hydrodynamically- 
assisted interactions, which could be much stronger than 
the pair vesicle interaction considered here. Thus, the ex- 
periments provide a lower bound for the interaction scale 
at dmin/R ~ 3.2 3.7, and so the corresponding upper 
bound for the volume fraction of non-interacting vesicle 
suspension is <f> — 0.08 -j- 0.13. 
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